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ABSTRACT 

The detection of upward propagating internal gravity waves at the base of the Sun's 
chromosphere has recently been reported by Straus et al., who postulated that these 
may efficiently couple to Alfven waves in magnetic regions. This may be important 
in transporting energy to higher levels. Here we explore the propagation, reflection 
and mode conversion of linear gravity waves in a VAL C atmosphere, and find that 
even weak magnetic fields usually reflect gravity waves back downward as slow mag- 
netoacoustic waves well before they reach the Alfven/acoustic equipartition height at 
which mode conversion might occur. However, for certain highly inclined magnetic 
field orientations in which the gravity waves manage to penetrate near or through 
the equipartition level, there can be substantial conversion to either or both upgo- 
ing Alfven and acoustic waves. Wave energy fluxes comparable to the chromospheric 
radiative losses are expected. 

Key words: Sun: atmospheric motions; Sun: magnetic fields. 



1 INTRODUCTION 

Using the Interferometric BIdimensional Spectrometer 
(IBIS) and the Echelle Spectrograph on the Dunn Solar Tele- 
scope (DST) of the Sacramento Peak National Solar Obser- 
vatory, and the Michelson Doppler Imager (MDI) on the 
Solar and Heliospheric Observatory (SOHO), Straus et al 

(2008) have identified upward propagating 1 gravity waves 
with frequencies between 0.7 mHz and 2.1 mHz in weak 
magnetic field regions of the solar atmosphere, and found 
that their energy flux is an order of magnitude larger than 
co-spatial acoustic waves, and comparable to the expected 
quiet-sun chromospheric losses of around 4.3 kWm~ 2 . How- 
ever, the gravity waves were found to be significantly sup- 
pressed in stronger field regions. Nevertheless, in light of re- 
cent identification of ubiquitous Alfven waves in the corona 
(De Pontieu et al. 2007; Tomczyk et al. 2007), they postulate 
that when these gravity waves enter magnetic regions they 
may efficiently couple to Alfven waves, perhaps contributing 
to the observed coronal wave flux. Furthermore, Jess et al. 

(2009) have directly identified torsional Alfven waves in Ha 
bright-point groups at frequencies as low as 1.4 mHz, which 
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1 The group velocity and hence energy flux is upward. As ex- 
pected of gravity waves, the phase velocity is downward. 



places them at least partially in the gravity wave regime if 
there is any coupling. 

In this paper, we explore the propagation, reflection, 
and mode conversion of atmospheric gravity waves of around 
1 mHz in frequency using both dispersion relations and nu- 
merical solution of the governing differential equations in 
simple atmospheric models with uniform inclined magnetic 
field. The imposed fields of 10 to 100 Gauss are weak in 
the photospheric context, but become dominant at greater 
heights as the plasma /3 (the ratio of plasma to magnetic 
pressure) falls below unity due to density stratification. We 
therefore have a situation where low frequency waves are 
essentially gravity waves at low altitudes, but become mag- 
netically dominated at higher levels in the chromosphere. 
The central question is: What happens to upward propa- 
gating gravity waves as they enter regions where magnetic 
forces become significant? Does the magnetic field help or 
hinder propagation through the chromosphere? The answer 
is: both, depending on magnetic field orientation. 



2 MODEL AND EQUATIONS 

We adopt the horizontally invariant VAL C Model of Ver- 
nazza et al. (1981), as adapted by Schmitz & Fleck (2003), 
up to z = 1.6 Mm. The interesting wave reflections and con- 
versions all happen below this height. An isothermal top is 
appended above 1.6 Mm. No transition region is included, 
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to avoid the complication of reflections outside our region of 
interest. 

The linear adiabatic wave equations for this scenario 
are set out in Cally & Goossens (2008), where exact series 
solutions are derived for the isothermal case. These are used 
here to specify top boundary conditions, which are always 
applied in the isothermal layers. 

Without loss of generality, the waves we consider are as- 
sumed to propagate in the x-z plane. The orientation of the 
magnetic field is then adjusted to explore various geome- 
tries. We distinguish the two-dimensional (2D) and three- 
dimensional (3D) cases: 

2D: The magnetic field lies in the same vertical (x-z) plane 
as the direction of wave propagation, and is inclined an angle 
9 from the vertical; 

3D: The vertical plane containing magnetic field lines 
makes an angle (f> 7^ to the x-z plane. 

In general, in cartesian coordinates (x,y,z), the magnetic 
field vector is B = B (sin 6 cos <j>, sin#sin<^, cos 6), where 
B = \B\. In light of the recent identification of ubiquitous 
horizontal magnetic field of up to 50-100 G in quiet inter- 
network regions using the Solar Optical Telescope/Spectro- 
Polarimeter (SOT/SP) on Hmode (Lites et al. 2008), par- 
ticular attention shall be given to highly inclined field. 

We use two mathematical tools: the magneto-acoustic- 
gravity dispersion relation (and the ray calculations which 
derive from it), and numerical integration of the 4 th order 
(2D) or 6 th order (3D) wave equations. 

2.1 Dispersion Relation 

Dispersion relations are widely used to describe oscillations 
in weakly inhomogeneous media, where the properties of the 
medium vary slowly on the length scale of a typical wave- 
length. This property allows us to effectively Fourier analyse 
in space (and time assuming a steady or slowly varying back- 
ground) to obtain a relationship between frequency lj and 
wavevector fe. 

The development of dispersion relations is intimately 
connected with determination of where the solutions of 
wave equations will be oscillatory and where non-oscillatory. 
This is a difficult enough task for gravitationally strati- 
fied unmagnetized atmospheres (Schmitz & Fleck 2003), 
let alone magneto-atmospheres. The difficulty ultimately re- 
volves around the definition of the acoustic cutoff frequency 
lo c which, depending on choice of dependent and indepen- 
dent variables, can take many distinct forms. There is simply 
no such thing as the acoustic cutoff frequency, and hence 
no firm and unambiguous way to precisely distinguish os- 
cillatory and non-oscillatory regions. The most commonly 
quoted expression is that of Deubner & Gough (1984), 

where c is the sound speed and H is the density scale 
height. Height z increases upward. However, as pointed out 
by Schmitz & Fleck, the appearance of the second deriva- 
tive of the density in this expression renders it largely im- 
practical for use with tabulated atmospheres in many cases. 
With the widely used Model S (Christensen-Dalsgaard et al. 
1996), it yields an enormous but very thin spike just below 
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Figure 1. Isothermal acoustic cutoff u) c i (thick full) and Brunt- 
Vaisala (dashed) frequencies as functions of height z for the VAL 
C model atmosphere. The Deubner & Gough cutoff frequency 
u) d g (thin full) is shown for comparison. 



the solar surface (see fig. 1 of Schunker and Cally 2006), 
which is inconsistent with the assumption of slow variation 
of coefficients for any solar waves of practical interest. 

Although Schmitz & Fleck (2003) derive more attrac- 
tive alternatives in the hydrodynamic case, these have not 
yet been extended to magnetohydrodynamics (MHD). We 
are therefore in uncertain territory applying dispersion rela- 
tions to solar atmospheric waves, especially in the few hun- 
dred kilometres below the surface. However, above the sur- 
face, the situation is less troublesome. All suggested formu- 
lations for uj c agree in the isothermal atmosphere case, where 
dH/dz = 0, returning the so-called isothermal acoustic cut- 
off frequency 



Indeed, for an isothermal atmosphere, there are exact so- 
lutions which confirm precisely this expression: see Lamb 
(1932, §309) for the hydrodynamic case and Cally (2001) 
with uniform magnetic field added. Consequently, in the 
lower solar atmosphere, where H varies (comparatively) lit- 
tle and slowly, we may adopt equation (2) with the expecta- 
tion of obtaining at least qualitatively correct results. Fur- 
thermore, for the gravity waves which interest us most here, 
it is the Brunt- Vaisala frequency N (about which there is 
no dispute) that is most relevant. 

Figure 1 plots the acoustic cutoff and buoyancy fre- 
quencies in the model atmosphere. Most significantly, the 
Brunt- Vaisala frequency in the VAL C atmosphere decreases 
steeply from around 6 mHz to approximately 1 mHz between 
0.6 Mm and 1.5 Mm, indicating a propensity to trap gravity 
waves in the low chromosphere. 

The 3D MHD dispersion function 

-T-, 2 2 2 7 2 1 / 2 2 , 2\ „ 

V = uj u> c a y k h + (lj — a fey ) X 

r 4 / 2 . 2\ 2 , 2 . 22,2,2 

[lj — (a + c )lj k + a c k k\\ 

+ c 2 N 2 k^ — (to 2 — a 2 k 2 )LJ 2 ] , (3) 

is derived in Appendix A. Here an exp[i(fc-;c — ojt)] depen- 
dence on position x and time t is assumed, k h — k x and 
k\\ are respectively the horizontal and field-aligned compo- 
nents of the wavevector, k — \k\ is the wavenumber, a is 
the Alfven speed, a y is the component of the Alfven veloc- 
ity in the y-direction (perpendicular to the vertical plane of 
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Figure 2. k z vs. z dispersion diagram for the non-magnetic case 
at 1 mHz with k x = 2 Mm -1 . The full curve represents T> = 
with the isothermal acoustic cutoff formula, and the dashed 
curve is for the Deubner & Gough formula. The lower branch 
represents a gravity wave with upward energy flux (as indicated 
by the arrow). Reflection occurs at z ~ 1.2 Mm for this frequency 
and horizontal wavenumbcr. 



propagation) and a z its vertical component, N is the Brunt- 
Vaisala (buoyancy) frequency given by 



iV 2 = 



9_ 
H 



(4) 



where g is the gravitational acceleration, and lo c = c/2H 
is the acoustic cutoff frequency. The dispersion function (3) 
reduces to that set out in equation (12) of Schunkcr and 
Cally (2006) in the 2D case = 0, with the addition of the 
decoupled Alfven factor lo 2 — a 2 fey . It also takes the usual 
magnetoacoustic form lo 4 — (a 2 + c 2 )uj 2 k 2 + a 2 c 2 k 2 k 2 in the 
absence of stratification, and the standard acoustic gravity 
wave form lo 2 (lo 2 — lo 2 ) + c 2 N 2 k 2 l without magnetic field. 
The first term L0 2 L0 2 ayk 2 provides the coupling between the 
magnetoacoustic and Alfven waves in three dimensions. 

The dispersion relation V = restricts allowable solu- 
tions in k-ui space. Now, in a horizontally and temporally 
invariant medium as is assumed here, lo and the horizontal 
component k x of k do not change with height, but the verti- 
cal component k z does. We may therefore fix lo and k x and 
plot dispersion curves in z-k z space. This gives a very infor- 
mative overview of the propagation properties of the various 
modes, and in particular where they reflect. 

Figure 2 displays a representative propagation diagram 
for the VAL C atmosphere in the absence of magnetic field. 
As there is little difference in the curves with u> c i and w D g, 
the former will be used exclusively throughout. The disper- 
sion relation reduces to 

jy 2 , .2 . .2 , .2 



kl 



k x + 



(5) 



in this case, clearly indicating vertical evanescence for N < 
lo < lo c for all k x . Even if lo < N < lo c , the wave is evanescent 
for small enough k x . This explains why the gravity wave 
dispersion curve in Fig. 2 is trapped below about 1.2 Mm in 
a "gravity cavity" . 

It is important to note that the upgoing branch in Fig. 2 



is the lower one, given the well-known property of gravity 
waves that the vertical components of the phase and group 
velocities have opposite signs (Lighthill 1978). Thus grav- 
ity waves with positive upward energy flux have downward 
phase velocity, i.e., k z < 0. 



2.1.1 Low 13 Asymptotics 

The asymptotic behaviour of wave modes in the o»c region 
is of crucial importance. A wide description of the behaviour 
of slow waves in this regime is given by Roberts (2006), in- 
cluding a discussion of applications to coronal loops. Mag- 
netically structured media [e.g., flux tubes) are beyond the 
scope of this paper, and will not be considered here. 

Analysis of the dispersion relation in the low-/3 regime 
reveals that there are two (possible) oscillatory classes of 
solution: 



Field-aligned acoustic (slow) waves lo 2 
os 2 9, i.e., 

yjLO 2 — L0l COS 2 9 



c 2 kf, + 



— ta,n6 cos <f>k x as a — > oo, (6) 
c cos 9 

provided lo > lo c cos 9. We term the reduction of the effective 
acoustic cutoff frequency by the factor cos# the "ramp ef- 
fect" . Equation (6) agrees perfectly with the exact /13 Frobc- 
nius eigenvalue in equation (26) of Cally & Goossens (2008). 
• The Alfven wave (which is also field-aligned of course) , 



2 2,2. 
lo ~ a fen + 



lo 2 — lo 2 cos 2 9 1 + tan 2 9 cos 2 < 



1 - 



lo 2 — lo 2 cos 2 6 1 + tan 2 9 cos 2 <j> 

— tan 9 cos <f>k x , (7) 



also as a — > 00. 



In both cases, the term k z = — tan 9 cos <f> k x alone would 
keep the phase of exp[i(fc x a; + k z z)] constant on field lines 
x — z tan 9 cos 0+constant, so it is the remaining square root 
terms which determine phase along field lines for these two 
field-aligned wave modes. In other words, — tan 9 cos <f> k x 
disappears if z is redefined as height along a fixed field line 
rather than at fixed x. Such a term is expected in the low-/3 
regime where the rigidity of the magnetic field lines domi- 
nates. It is purely geometric. 

The remaining fast magneto-acoustic wave is always 
evanescent for fixed lo and k x and large enough a»c since 
its dispersion relation is fc 2 w lo 2 /a 2 — k 2 < . 



2.2 Numerical Integration 

Although convenient and informative, the dispersion rela- 
tion is only an approximate tool. We corroborate our results 
by numerically integrating the linear adiabatic wave equa- 
tions across < z < 1.6 Mm. The integrated equations do 
not depend upon lo c explicitly, and so are not subject to the 
uncertainties associated with that term. The top boundary 
conditions are taken from the exact convergent series so- 
lutions derived by Cally & Goossens (2008) for isothermal 
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slabs, where a 3> c, and which represent (i) an outgoing 
or evanescent slow (acoustic) wave (depending on whether 
u> > u) c cos# or vice versa; (ii) an evanescent fast (magnetic) 
wave; and (iii) an outgoing Alfven wave (in the 3D case). 
There are no incoming or exponentially growing solutions 
allowed at the top. 

At the bottom (z ^ 0) we generalize the radiation 
boundary condition (8) of Cally & Goossens (2008) to suit 
gravity waves as well as acoustic waves. The aim is to disal- 
low incoming slow (magnetic) and Alfven waves at the base, 
where a <C c, but allow incoming gravity or acoustic waves, 
and outgoing waves of all varieties. We do this by imposing 

9 -Wo ■ fd 



dt 



- a-V 



and 



(£—v)^ = 0,(8) 



where £ = (£, 77, £) is the plasma displacement vector 
(though the r\ equation is redundant in 2D where 77 decou- 
ples), and A is the gravito-acoustic annihilator operator de- 
fined by 



p 



-1/4 



N 2 £ . , N 2 C 
g u 2 



(9) 



The square bracket term in equation (9) is identically zero 
in the absence of magnetic field. Therefore, in the region 
a « c where the gravito-acoustic and magnetic waves are 
fully decoupled, it strongly suppresses the gravito-acoustic 
waves, leaving only the slow (magnetic) waves. The factor 
p _1//4 in (9) removes an amplitude factor in the remaining 
sinusoidal slow waves. The d/dt — a-'V operator then selects 
the upcoming magnetic wave, which is set to zero with this 
condition. The incoming Alfven wave is similarly suppressed 
by the second equation in (8). 2 

The differential equations being of 4 th (6 th ) order in 2D 
(3D), and with 2 (3) boundary conditions applied at the top 
and 1 (2) at the bottom, an arbitrary normalization condi- 
tion is all that remains to be applied. The gravity waves are 
assumed generated at z — (by granulation), so it is not 
necessary to impose an evanescence condition on them be- 
low this level, which would result in an eigenvalue problem, 
rather than the driven wave scenario we envisage here. 

As a test of the bottom boundary condition, we apply 
it at several depths (z = 0, — 1, —2 Mm), and find that our 
solution typically varies only in the third significant figure. 
This indicates that equations (8) represent a good approx- 
imation to the required radiation condition. One could do 
better by developing an asymptotic series solution (as in 







Figure 3. Ray plots for rays of frequency 1 mHz and horizontal 
wavenumber k x = 2 Mm' 1 with (from top to bottom): (i) no 
magnetic field; (ii) vertical 10 G field; (iii) 10 G field inclined 80° 
from the vertical (2D); (iv) 10 G field inclined -80° from the 
vertical; (v) 10 G field inclined 80° from the vertical and oriented 
50° away from the x-z plane (3D); and (vi) x-y projection of the 
previous case. Field line orientation is indicated by the faint grey 
backgrounds (projected in case (v)). The dots on the ray curves 
indicate 1 minute group travel intervals. The a = c equipartition 
level is indicated by a horizontal line where it falls within the 
depicted domain. In all cases the rays are propagating rightward. 



2 The p^ 1 / 4 factors may be unfamiliar, and so warrant some ex- 
planation. The Alfven wave equation takes the form d 2 r)/dt 2 = 
a 2 d 2 r)/ds 2 , where s = zsec8 is distance along the field line and 



a oc p 



-1/2 



is the Alfven speed. Assuming harmonic time de- 
pendence exp(— hot) reduces this to d 2 rj/dz 2 = Q(z)r), where 
Q = — u> 2 /(a 2 cos 2 9) oc p. The WKBJ amplitude factor for such 
an equation is then |Q| — <x p~ x / 4 (Bender & Orszag 1978). 
Thus the p 1 / 4 factor in (8) expunges this amplitude dependence 
to leading order. The slow wave in a « c takes essentially the 
same form, though the polarization is in the x rather than y di- 
rection, and £ may be used as the dependent variable instead of 
r\. The dominant term in the square bracket in A£ is obviously 
d£/dz oc k z £. But k z ~ u)/acos6 oc p 1 / 2 and the amplitude of 
£ varies as p -1 / 4 as before, showing that d£,/dz has amplitude 
factor p 1 / 4 to leading order. The p" 1 / 4 factor in (9) cancels this. 
These arguments are confirmed numerically. 



Crouch & Cally 2005) and applying it at sufficient depth, 
but this is both messy and unnecessary, it would be specific 
to the particular model atmosphere in which it were applied, 
and it could not easily be done for a tabulated model. 



3 RESULTS 

For simplicity, we focus on 1 mHz waves with horizontal 
wavenumber k x — 2 Mm -1 , placing them well within the 
high flux regime identified observationally by Straus et al 
(2008). Two magnetic field strengths, B = 10 G and B = 
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100 G will suffice to illustrate how the gravity waves are 
influenced by typical chromospheric fields. 



3.1 Dispersion and Ray Diagrams 

Ray plots 3 for four cases are collected in Fig. 3. They illus- 
trate, respectively, 

(i) In the absence of magnetic field, the gravity wave re- 
flects symmetrically at around z = 1.2 Mm, as expected 
from Fig. 2. 

(ii) In vertical field, the wave reflects at a much lower 
height, and as a slow magneto-acoustic wave (we can see 
that it is "slow" by the very close spacing of the 1 min 
dots) . 

(iii) At high field inclination in 2D, the gravity wave con- 
verts to an acoustic wave which then follows the field lines. 
The regular grid spacing above a = c indicates that it is 
acoustic (the atmosphere is nearly isothermal). 

(iv) At high inclination in the other direction (negative 
9) the ray again reflects as a slow wave following the field 
lines. 

(v) At high field inclination in 3D, with <j> = 50° here, the 
mode conversion is to the Alfven wave, as is apparent from 
the rapidly increasing 1-min dot separation with height. The 
depicted ray path is a projection onto the x-z plane. In real- 
ity, this last ray turns sharply out of the x-z plane to follow 
the field lines, as illustrated in the bottom panel. 

Clearly, the determining characteristic here for the different 
behaviours is the "attack angle" (Schunker and Cally 2006) 
the ray makes to the magnetic field lines. At small attack 
angle the ray may avoid reflection. 

These are "classical" ray calculations (Weinberg 1962), 
as distinct from the generalized ray formalism of Schunker 
and Cally (2006) (recently tested against exact solutions by 
Hansen & Cally 2009). Therefore, they do not account for 
partial transmission/conversion. Instead the rays simply fol- 
low the connectivity implied by the dispersion curves. Never- 
theless, these are often (but not always) the dominant paths. 
Partial mode conversions are addressed in Section 3.2 using 
full numerical solution. 

Fig. 4 displays the dispersion diagrams corresponding 
to the ray paths depicted in Fig. 3, parts (ii)-(v) respec- 
tively. The non-magnetic case (part (i)) has already been 
covered in Fig. 2. Points to note include (with numbering 
corresponding to that of the previous list) 

(II) The dispersion curve turns over at a very low height 
in vertical magnetic field, even before a 2 /c 2 reaches the 1% 
level, in complete agreement with the corresponding ray 
path. The "return path" clearly corresponds to very slow 
propagation, as \k z \ becomes very large (recall that the 
phase speed is u)/\k\ and ui and k x are fixed). 

(III) With highly inclined 2D field the upgoing grav- 
ity branch penetrates the a = c equipartition level and 
quickly settles onto the low-/3 asymptotic acoustic curve 

3 The ray equations are derived from the dispersion function, 
which takes the role of a Hamiltonian. See for example Cally 
(2009a), equation (14). They are solved as a set of ordinary dif- 
ferential equations along their arc lengths, starting at x = 0, 
z = 0.1, just inside the gravity cavity. 
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Figure 5. Acoustic (full) and magnetic (dashed) wave energy flux 
(kWm~ 2 ) as functions of magnetic field inclination 9 and three 
orientations (0°, 30°, and 60° respectively) for B = 10 G and 1 
mHz waves with k x = 2 Mm -1 . In all cases, the vertical velocity 
w at z = is normalized to 1 kms -1 . Despite appearances, the 
rapid oscillations at high inclination are smoothly resolved, and 
arc not numerical artefacts. 



(blue dashed lines). Again, this is in total accord with the 
ray figure. 

(IV) When the highly inclined (80°) field is oppositely 
oriented (—80°) and the attack angle is large, there is no 
such connection to the low-/3 acoustic waves, and the gravity 
wave reflects once again as a slow magneto-acoustic wave, 
as indicated by the arrows in the lower left panel of Fig. 4. 
Note that this panel is simply the reflection about k z = of 
the top right panel. 

(V) In the 3D case the Alfven wave now couples to the 
acoustic-gravity waves (the Alfven loci have been suppressed 
in the first three panels, as they were inaccessible in 2D) . In 
the case depicted, the upgoing gravity wave connects directly 
to the Alfven wave, as indicated by the red dotted curve 
representing the Alfven asymptotics of equation (7). Once 
again, the agreement with the ray figure is apparent. 



3.2 Transmitted Fluxes 

The wave-energy flux carried by linear MHD waves may be 
broken into gravito-acoustic and magnetic parts, F = f ac + 
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Figure 4. Dispersion diagrams for waves of frequency 1 mHz and horizontal wavenumber k x =2 Mm . The vertical lines indicate the 
value of the ratio a 2 /c 2 at various heights: thick solid line - a 2 /c 2 = 1 (the equipartition level); dashed line - a 2 /c 2 = 0.1; dotted line - 
a 2 /c 2 = 0.01. Top left: vertical 10 G field - the presence of even a weak vertical magnetic field causes the up-going gravity wave (bottom 
branch) to reflect back downwards as a slow magneto- acoustic wave, and it does so even before a 2 /c 2 reaches 1% of the equipartition 
value; Top right: 10 G field inclined 80° from the vertical (2D) - when the field is highly inclined and the attack angle is small, the 
up-going gravity wave penetrates the equipartition level and connects to the up-going field-aligned acoustic wave (given asymptotically 
by equation (6) which is indicated here by the blue dashed lines); Bottom left: 10 G field inclined —80° from the vertical — when the 
attack angle is large, the up-going gravity wave is reflected down as a slow magneto-acoustic wave (as indicated by the arrows); Bottom 
right 10 G field inclined 80° from the vertical and oriented 50° away from the x-z plane (3D) - in this case the up-going gravity wave 
connects to the Alfven solution (asymptotic equation (7), the red dotted curve). These results are in complete agreement with the ray 
behaviours of Fig. 3. 



J 1 mag, where F ac = Re[piv*] and F mag = -Re[J5i x2?*], 
pi is the Eulerian gas pressure perturbation, v is the plasma 
velocity, and Ei and B 1 are the electric and magnetic field 
perturbations. Once the numerical solution is obtained, the 
vertical components F ac and Fmag may be calculated in the 
isothermal layer at the top (see Cally 2009b, for details). In 
3D, the flux of most interest is Fmag, the Alfven flux escap- 
ing from the top (identified as such since the fast wave is 
evanescent in that regime). In all cases, the vertical velocity 
w at z = is normalized to 1 kms" 1 . Since the flux scales 
as the square of the velocity perturbations, these results are 
easily adjusted for different driving velocities. 



Fig. 5 for B = 10 G shows substantial acoustic flux 
in the 2D case cj> — 0° , once the ramp effect has turned 
on (cosd < uj/uj c ). In practice, this requires highly inclined 
field, characteristic of canopy. As the field is rotated out of 
the x-z plane, coupling to the Alfven wave takes over, and 
the acoustic flux diminishes. However, again it is only signif- 
icant at high field inclination. Very similar results are found 
at 100 G (Fig. 6). It is most interesting that the resultant 
fluxes, even with a modest base driving velocity of w = 1 
kms" 1 , are comparable to the estimated chromospheric ra- 
diative losses, and are therefore of genuine significance. 

It is important to not give too much weight to the wave 
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Figure 6. As for Fig. 5, but with B = 100 G. 




connectivities implied by the eikonal approximation through 
the dispersion and ray diagrams. As is familiar from quan- 
tum mechanics, energy can jump gaps between (closely) 
neighbouring dispersion curves. This is well illustrated by 
Fig. 7. The second panel of Fig. 5 shows that both acoustic 
and Alfvenic flux exits the top in this 3D case, though with 
the magnetic flux dominating. Comparing with the disper- 
sion curve, the eikonal connectivity is to the Alfven branch, 
but the gap (avoided crossing) to the acoustic branch is nar- 
row. This is clearly indicative of partial tunnelling across the 
gap. 

Fig. 8 fixes 6 = 85° and rotates </> from 0° to 170°. Once 
again, we see that the acoustic flux diminishes as the field 
orientation moves away from the 2D configuration, but that 
it is largely compensated by an increase in Alfven flux, at 
least out to about 70°, after which both fluxes quickly be- 
come insignificant. This is reminiscent of Fig. 2 in Cally & 
Goossens (2008), though shifted here to much higher 9 be- 
cause of the lower frequency. Fig. 9 illustrates this further 
using dispersion diagrams. In 2D the Alfven branch is in- 
accessible, and so carries no flux. However, as the magnetic 
field is rotated out of the x-z plane, the upgoing gravity 
wave branch instantly connects to the Alfven branch at large 
z. Most energy continues to flow along the low-/3 acoustic 
branch though by jumping the (small) gap that has opened 
up as an avoided crossing (primarily because the polariza- 
tions are still almost orthogonal at small <j>). At larger <f> (60° 
and 74° are illustrated here) the gap is very large and so the 
energy flows overwhelmingly along the dispersion curve onto 
the Alfven branch. However, by cj> = 80° this upgoing Alfven 
branch has been cut off (due to the term under the square 
root in (7) becoming negative) and the connectivity is now 
to a downgoing Alfven wave, explaining the precipitous drop 
in Alfven flux with increasing <f> seen in Fig. 8. 

Fig. 10 illustrates how wave energy flux may shift from 
acoustic to magnetic near the a — c equipartition level in 3D. 
Also note the negative magnetic flux at the bottom, indicat- 
ing some significant reflected Alfven and/or slow magneto- 
acoustic flux there. Typically, the conversion to Alfven waves 
occurs over a broader and slightly higher region than is 
characteristic of fast/slow conversion, and that is apparent 
here, with a steady conversion process taking place between 
roughly 1 Mm and 1.5 Mm. The fact that the dispersion lo- 
cus which penetrates a — c in Fig. 7 essentially sits right on 
top of the asymptotic Alfven curve (red dotted) shows just 
how well the acoustic-gravity and Alfven phase velocities 
match over an extended interval, which is very favourable 
for mode conversion (Cally 2005). 



Figure 7. Dispersion diagram for a wave with frequency 1 mHz 
and horizontal wavenumber k x =2 Mm -1 , subject to a 10 G field 
with 9 = 83°, <t> = 30°. The gravity wave dispersion curve con- 
nects to the Alfven wave above a = c. However the field guided 
acoustic wave (asymptotic to the upper blue dashed line) passes 
close to this gravity /Alfven locus near the equipartition level, and 
energy can tunnel across the gap (c.f. the middle panel of Fig. 5). 



Fig. 11 briefly addresses how the escaping wave energy 
fluxes at the top vary with frequency. The two cases shown 
are for 0.7 mHz and 2.1 mHz, but otherwise address the same 
case as in the middle panel of Fig. 5. The results are qual- 
itatively similar, except that the flux "turns on" at greater 
inclinations for lower frequency, and lesser inclinations for 
higher frequency, as expected from the ramp effect, which 
allows acoustic waves to propagate once ui > LJcCos^. It is 
also notable that the magnitude of the Alfvenic flux is con- 
siderably increased at higher frequency. 
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1 mHz, k,=2 Mm"', 8=85°, B=100 G 
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Figure 8. Acoustic (full) and magnetic (dashed) wave energy flux (kWm~ 2 ) as functions of magnetic field orientation <f> with inclination 
8 = 85°. Left panel: B = 10 G; right panel: B = 100 G. As in figures 5 and 6, the frequency is 1 mHz and the horizontal wavenumber is 
k x = 2 Mm -1 . The vertical velocity w at z = is normalized to 1 kms" 1 . 




Figure 9. Dispersion diagrams illustrating the role of field orientation <j> in the waxing and waning of Alfven energy loss. All cases 
correspond to field inclination = 85°, field strength 10 G, frequency 1 mHz, and horizontal wavenumber 2 Mm -1 . Top left: <j> = 0°, 
with the Alfven branch greyed out to denote that it is decoupled from the gravito-magneto-acoustic branches; Top centre: <fi = 5°; Top 
right: <f> = 60°; Bottom left: <j> = 74°; and Bottom centre: <j> = 80°. 



4 CONCLUSIONS 

The main conclusions we draw from our analyses are: 

1. Even very weak magnetic fields very effectively re- 
flect gravity waves back downward as slow magneto-acoustic 
waves. This typically happens well below the a = c equipar- 
tition level. 

2. However, at very large magnetic field inclinations, typ- 



ically around 80° or more depending on frequency (see 
Fig. 11), substantial mode conversion from gravity waves 
to either field-guided acoustic waves (for small (j>) or Alfven 
waves (20° < (j> < 70°) occurs, and these waves continue to 
propagate upward along the field lines. The amount of en- 
ergy they carry is potentially significant for the upper chro- 
mosphere. 

3. Wave energy fluxes reaching the top of our model are 
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very sensitive to magnetic field direction, but quite insensi- 
tive to magnetic field strength, at least in the range 10 G - 
100 G. 

4. The dispersion diagrams give a simple, easy and quite 
accurate picture of the behaviour of gravity waves in a 
magneto-atmosphere, though with the caveat that tun- 
nelling can sometimes occur between branches. 



In simple terms, we conclude that atmospheric grav- 
ity waves are very effectively suppressed by even very weak 
magnetic field, unless that field is highly inclined and the 
attack angle fine, in which case it opens a window to the 
upper atmosphere that allows the gravity waves to propa- 
gate through in a different guise. This is closely related to 
the "magnetic portals" of Jefferies et al. (2006) which rely on 
the ramp effect to allow acoustic waves below the acoustic 
cutoff frequency to still propagate upward in low-/? inclined 
magnetic field, but in the case of low-frequency gravity waves 
which are already propagating, it gives them the opportu- 
nity to convert to propagating acoustic waves around the 
a = c level. Because of their low frequency though, very 
substantial inclination is required to open these windows, 
but no more than is characteristic of chromospheric canopy. 

There are consequences for recent and future observa- 
tions of solar atmospheric oscillations. The surprising extent 
to which even very weak vertical or moderately inclined 
magnetic field inhibits gravity waves by causing them to 
quickly reflect as slow magneto-acoustic waves perhaps ex- 
plains "significantly suppressed atmospheric gravity waves 
at locations of magnetic flux" found by Straus et al (2008). 
The ubiquity of near-horizontal field in the low solar atmo- 
sphere discovered recently with Hinode (Lites et al. 2008) 
however raises the possibility that low frequency gravity 
waves may efficiently couple to Alfven waves that continue 
to propagate vertically into the corona, contributing to the 
vast sea of waving field lines now known to exist there (De 
Pontieu et al. 2007; Tomczyk et al. 2007). Although Tom- 
czyk et al. detect a peak in velocity power of these Alfven 
oscillations at around 3.5 mHz, (presumably driven by the 
Sun's internal normal modes), the power spectrum contin- 
ues to rise with decreasing frequency till at least 1 mHz, 
well inside the gravity wave regime at photospheric level. It 
is tempting to postulate that gravity waves may be the vec- 
tor of this wave energy at low levels and that it may convert 
to Alfven waves around the acoustic/ Alfvenic equipartition 
level in highly inclined field regions. Further observational 
work, ideally at multiple heights, is warranted to more fully 
explore these possibilities. 

It should be emphasised though that our analysis is 
entirely linear. Acoustic waves are likely to shock before 
reaching the upper chromosphere. However, Alfven waves 
do not suffer this fate. Our models are also adiabatic, which 
is not a good representation of the chromosphere, though the 
detections of Straus et al (2008) suggest that atmospheric 
radiative losses do not completely suppress gravity waves, 
at least at the photospheric altitudes sampled by IBIS and 
MDI. The adiabatic assumption will be relaxed in future 
work. 



Sa/\/\A 
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Figure 10. Acoustic (full) and magnetic (dashed) wave energy 
fluxes as functions of height z (Mm) for B = 10 G, 6 = 83°, 
<j> = 30° , corresponding to Fig. 7. Once again, the frequency is 1 
mHz and the horizontal wavenumber is 2 Mm -1 . The vertical line 
indicates the position of the Alfvcn/acoustic equipartition level 
a = c. An extended section of the isothermal 
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Figure 11. Acoustic (full) and magnetic (dashed) wave energy 
flux (kWm~ 2 ) as functions of magnetic field inclination 9 with 
(f> = 30° for B = 10 G, k x = 0.2 Mm -1 , and two different fre- 
quencies. Top panel: 0.7 mHz; Bottom panel: 2.1 mHz. The middle 
panel of Fig. 5 is intermediate between these two cases. Note the 
different 0-scalcs on the two graphs. 
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APPENDIX A: DISPERSION RELATION 

As explained in Schunker and Cally (2006), a convenient 
starting point for deriving the dispersion relation in MHD is 
the Lagrangian density associated with the linearized MHD 
equations (Goedbloed & Poedts 2004). A little manipulation 
using the equilibrium condition Vp = pg + jxB yields the 
elegant Hermitian form 



where £ is the plasma displacement vector and £ its time 
derivative, p the density, c the sound speed, B the mag- 
netic field, j = V X B the current density, g — —ge z the 
gravitational acceleration, and b — Vx (^xB) the mag- 
netic field perturbation. The magnetic permeability p, has 
been scaled to unity for simplicity. Variation of the action 
/ C dV dt with respect to £ recovers the linear wave equa- 
tions. The last term in Equation (Al) vanishes for a force 
free magnetic field, and the last two terms for a potential 
field. For the case at hand, a uniform field, both terms are 
dropped. We choose to adopt X = p 1 ^ 2 ^ as the dependent 
variable in the non-magnetic terms, but to remain with £ in 
the magnetic term Vx(^xB). 

Following the standard eikonal method (Weinberg 
1962), we may then identify VX = ifeX, and in particu- 
lar V-X = ifc-X in the non-magnetic terms. On the other 
hand, Vx(^xB) = ifcx (Xx a) is adopted in the magnetic 
energy term. Assuming the background state is stationary, 
dX/dt ee -iwX. Then 



£=^ 2 |X| 2 



1 2 

2 C 



fe-X 



iZ 
2H 



-|fcx(Xxa)| 2 



+ ^|(Z*fc-X-Zfc-X*) 

= ^ 2 -a 2 fc 2 )|X| 2 -i(a 2 + c 2 )|fc.X| 2 

+ i (m-X (fc-Xr + (m-X)* fc-X) - ^lj 2 c \Z\ 2 , 

(A2) 

where Z is the vertical component of X, e z is the up- 
ward vertical unit vector, a = p~ 1,>2 B = a(sin 9 cos <j> e x + 
sm8sm<f>e y + cos0e z ) is the Alfven velocity (with a = a 
the Alfven speed), and k\\ = a-k = k cos a is the field-aligned 
component of k. Here, we have found it convenient to define 
the complex vector 



m = (a-k) a — i ( g 



2H 



(A3) 



From the Lagrangian density, we may read off the Her- 
mitian dispersion tensor D defined by C — |X H DX in ma- 
trix notation, and consequently its determinant, the disper- 
sion function, which simplifies to 

V = dct D = ui 2 LU 2 a 2 k 2 + (to 2 — a 2 k 2 ) x 

r4 / 2 , 2x2,2. 22,2,2 

|o> — (a + c )ui k + a c k fcy 

+ c 2 N 2 kl - (lu 2 - a 2 fc 2 K 2 ] , (A4) 

where a z is the vertical component of the Alfven velocity 
and a y is the component perpendicular to the plane con- 
taining k and g. The Brunt- Vaisala frequency iV is defined 
by iV 2 = g/H — g 2 /c 2 , lu c = c/2if is the isothermal acoustic 
cutoff frequency, and kh is the horizontal component of the 
wavevector. 
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